Massive Dirac fermions in single-layer graphene 
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Motivated by the results of recent photoemission and tunneling studies, we discuss potential 
many-body sources of a finite gap in the Dirac fermion spectrum of graphene. Specifically, we focus 
on the putative Peierls- and Cooper-like pairing instabilities which can be driven by sufficiently 
strong Coulomb and electron-phonon interactions, respectively. Our results compare favorably with 
the available experimental and Monte Carlo data. 



Introduction 

The recent advent of graphene has brought about a 
plethora of theoretical studies of the emergent behav- 
iors and new phases of electronic matter in degenerate 
semimetalic systems. 

Amongst other spectacular experimental findings, re- 
cent photoemission^ and tunneling^ data obtained, corre- 
spondingly, in epitaxially grown and suspended graphene 
demonstrated the presence of a sizable gap in the (othe- 
wise, nearly linear) spectrum of the Dirac-like quasiparti- 
cles which dominate the transport and optical properties 
of this novel material. Moreover, comparable gaps were 
observed in bi- and tri-layers. 

Thus far, the proposed explanations of these findings 
invoked the effect of a commensurate static corrugation 
(surface reconstruction, "frozen phonon") and/or break- 
ing of the symmetry between two sublattices (hereafter, 
A and B) of the graphene sheet due to a nearby sub- 
strate, which mechanisms can result in a hybridization 
of the electronic states from the vicinity of the two Dirac 
points (hereafter, L and i?)^'^'^. 

Alternatively, the gap can also be generated by the 
Coulomb or phonon-mediated interactions, thus resulting 
from an incipient instability of the nodal fermion system 
itself. 

In particular, the pertinent Coulomb interaction- 
driven Peierls-likc instabilities can be divided onto two 
main groups. The first one is represented by various 
graphene-specific modifications of the standard scenario 
of excitonic insulating behavior developing in fermion 
systems with a finite Fermi surface^. 

Specifically, such proposals were discussed in the cases 
of a parallel field^i^ and an electrically biased bi-layer^, 
under either of which conditions one has matching pairs 
of electron and hole Dirac gases with finite (and equal - 
in the case of undoped graphene) chemical potentials, as 
measured from the Dirac point. 

In these BCS-type scenarios, an opening of the exci- 
tonic gap at arbitrarily weak couplings is facilitated by 
a finite density of states (DOS) at the Fermi surface, the 
gap being proportional to the effective Fermi energy of 
either of the two mutually attracting components. In par- 
ticular, the authors of Ref.— found the excitonic gap to 
scale (almost) linearly with the Zeeman splitting, A ~ _B. 

In contrast to these scenarios, the possibility of exci- 
tonic pairing between the Dirac fermions with a vanishing 
DOS was first discussed in the context of multi-layered 



HOPG graphite, and later graphene, for both, the physi- 
cal long-range Coulomb interactional^ as well as its short- 
range surrogates'^. By and large, the justification for 
using the latter, thus far, has been their computational 
convenience, as compared to the genuine Coulomb case. 
The early analysis carried out in Refs^i^ predicts the 
existence of a finite threshold for the strength of the 
Coulomb coupling quantified in terms of the analog of 
the "fine structure constant" 
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where e is the dielectric constant of the surrounding 
medium (in the case of a graphene monolayer sandwiched 
between two different media it equals the average of the 
corresponding dielectric constants). 

Moreover, Refs*?'^ also discussed a concomitant phe- 
nomenon of 'magnetic catalysis', according to which 
a perpendicular magnetic field can open up a gap ^ 
y/B — Bq, where Bq is proportional to the electron den- 
sity of doping relative to the Dirac point, at arbitrar- 
ily small values of g^. Unlike the mechanism proposed 
in Ref.^, the gap-opening effect of the magnetic field 
appears to be of a purely orbital nature. Also, con- 
trary to the unproven experimental status of the pro- 
posal of Refi^, the field-induced gap predicted in Refsi^i^ 
may have already been observed in both, HOPG^i and 
graphene^^ (the over an order of magnitude disparity be- 
tween the measured critical fields Bq might stem from a 
comparable difference in the doping densities of the two 
systems). 

However, the feasibility of the pairing mechanism of 
Refs.^'^ was later questioned in^ii^ whose authors argued 
that it would likely be completely hindered by a down- 
ward renormalization of the effective Coulomb coupling 
(1) due to an increase of the renormalized Fermi velocity 
with decreasing momentuirii* . 

In what follows, we refine the analysis of Refs.^ and 
extend it to the cases of Coulomb- and phonon-mediated 
Peierls- and Cooper-pairings, respectively, while account- 
ing for the effects of the running Fermi velocity. 

Pairing instabilities and mass terms 

With an eye on a unified description of both, Peierls- 
and Cooper- types of pairings, we describe the Dirac 
fermion states as 16-component vectors 
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where the first, second, and third indices in each of the 
electron- and hole-like 8-component parts of this vector 
pertain to the sublattice (C = A,B), valley (n = L,R), 
and spin {a =t, i) degrees of freedom, respectively. 

The set of operators acting in the corresponding 
Hilbert space is spanned by the basis of 4'' matrices 
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The four consecutive factors in the above tensor product 
act in the Nambu (particle-hole), sublattice, valley, and 
spin subspaces, each of the four indices taking one of the 
four possible values (n, a,i,a = 0, 1, 2, 3). 

With the magnetic field and chemical potential in- 
cluded, the Dirac fermion Green function reads 

G(e, p) = [e - P3 ® K||p - M + sB) + S(e, p)]-^ (4) 

where v ~ lO^m/s is the velocity of the Dirac-like nodal 
excitations. The part of the self-energy S^ = S commut- 
ing with the Dirac Hamiltonian Hq = ucr||P {[Hq, E] = 0) 
accounts for the renormalization of the Dirac fermions' 
dispersion via the effective Fermi velocity v{p), while the 
anticommuting one "Ea — A {{Hq,A} — 0) incorporates 
a (possibly, momentum-dependent) fermion mass func- 
tion A(p) = TrAVie. 

Together, they give rise to the massive Dirac fermion 
spectrum 
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The list of putative order parameters includes various 
Peierls- and Cooper-like fermion bilinears 
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Among these bilinears are those which can be associ- 
ated with 4 different types of spin-singlet Dirac-type 
mass terms vlftp^vp^ classified according to their sym- 
metry properties under the parity (P— ): A{B) 
B{A);L{R) ^ R{L) and time reversal (T-): A{B) 
A{B);L{R) ^^ R{L) transformations: 

ll^'^a^i (g) Ti (g) SoV' = A^^La^Ra + B^^Aro, + h.C. 

il^'^a-i (g ra So'/' = Y^ (AL^^^ - s|„B»a) 

i=L,R 

ip'' (T3 (» To (» soip = ^ sgni{Al^Aia - BJ^Bia) 

i=L,R 



i^'^Cry, (g) T2 (g) Soil) = iA^^BRa - iBl^ARa + h. 
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where we choose sgnL{R) = ±1 (hereafter, all the sum- 
mations over the spin indices are implicit). The above op- 
erators (in the listed order) are P, T-even; P-odd, T-even; 



P-even, T-odd; and P, T-odd, respectively. By replacing 
So with S3 one can also construct spin-triplet counter- 
parts of Eqs.(7). 

Notably, the above order parameters include both, 
intra- , < A(P)]^ j^A{B)l,r >, and inter-node, < 

^(^)I r^{-^)r,l >, excitonic pairings with the total 

momenta and Kl — Kr — 2Kl = ~Kl, respectively. It 
was shown in Ref.'^ that the former mass terms manifest 
themselves as a charge density wave (CDW) representing 
the excess/deficit of the electron density which alternates 
between the two sublattices. In contrast, the latter terms 
correspond to the Kekule dimerization pattern which re- 
sults in the tripling of the unit cell and the ensuing equiv- 
alence of the two Dirac points in the reduced Brillouin 
zone. 

For comparison, the ordinary and staggered chemical 
potential/magnetic field terms correspond to the bilin- 
ears 

Y^S'-^ihsgnaKAlA^ + BlB^p), 

i 

J^r'^il, sgnaUlAp bIb,^) (8) 

i 

which only cause energy level shifts, but no gaps. 

Likewise, an attractive interaction can potentially gen- 
erate the following Major ana- type mass terms 

V'O'O <g n <g SoV' = X! {^^aA^I3 + BiaB^i3)s2'^ 
i=L,R 

V'o-o <g T2 (g soV' = X! sgni{AiaAip + BiaBif3)s2'^ 

i=L.,R 
i>ao g) To g) S2V' = (ALaBRa - BLaARa) 
^Oq (g Ta (g SoV' = {AhaBRa ~ BLaARa)s2^ (9) 

amongst which there are both, inter- node (uniform), 
< Alm.Br,l >, and intra-node, < A{B)l,rA{B)lm. >, 
Cooper pairings, the latter corresponding to a non- 
uniform (LOFF-type) state. 

Which of the above orderings (or a combination 
thereof) would have a higher chance to occur, depends 
on the parameter values of the overall effective electron 
interaction which is comprised of the direct Coulomb re- 
pulsion and the effective phonon-mediated attraction. 

Electron interactions and excitonic pairing 

At small transferred momenta the effective interaction 
in undoped (/i — 0) graphene immersed in a medium 
with e ~ 1, is dominated by the Coulomb repulsion which 
brings about a putative Peierls-like(excitonic) pairing in- 
stability. Contrary to the model problem of A'' ^ 1 Dirac 
fermion species studied in Ref.^ by means of the 1/A^- 
expansion, that of graphene features no small parameter. 

By analogy with the earlier analyses^ of a fully rela- 
tivistically invariant counterpart, dubbed chiral symme- 
try breaking, of the excitonic instability the latter can 
studied in the framework of the Schwinger-Dyson equa- 
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(10) 
written in terms of the dressed fermion Green, interac- 
tion, and vertex (F) functions. 

In the case of an ultra-violet divergent momentum inte- 
gral the very derivation of Eq.(lO) becomes problematic 
because of the lack of a proper account of non-universal 
local terms in the Hamiltonian. While other works— re- 
sorted to a phenomenological approach where such terms 
would be treated as free parameters, we focus on the sit- 
uation where the gap equation receives a dominant con- 
tribution from the momenta A/v ^ q <C A, thereby jus- 
tifying the derivation of Eq.(lO). Incidentally, such a be- 
havior does occur in those situations where the fermion 
DOS vanishes and/or the interaction function remains 
sufficiently singular at small momenta. 

With the fermion polarization 
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included, the effective interaction in Eq.(lO) reads 
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While the vertex function in Eq.(lO) remains largely un- 
affected by the interactions, F « 1—, the dressed fermion 
Green function has to incorporate, in addition to the 
mass term, the aforementioned renormalization of the 
Fermi velocity. The latter (and, concomitantly, the effec- 
tive momentum-dependent coupling g(j))) were approxi- 
mated in Refi^i^ by a power-law 
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within the range of momenta A/w < p < A, where A is a 
momentum cut-off of order the inverse lattice spacing. 

Computed in the l/A'^-expansion, the exponent as- 
sumes the value^i^ 
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which remains relatively small down to the physically rel- 
evant number TV = 4 of two-component (pseudo)spinors 
with different values of the physical spin and valley in- 
dices. At still lower energies and/or go ^ 1 the renor- 
malized Fermi velocity shows an even slower (only loga- 
rithmic) momentum dependences,. 

Upon neglecting the unimportant vertex corrections, 
switching to imaginary frequences, performing the fre- 
quency integration, and putting s = vp in Eq.(lO), one 
arrives at the equation for the momentum-dependent 
mass function defined as the anticommuting part of the 



self-energy A{p) = T,a{ivp,p) evaluated on the mass- 
shell 



Mp) = Yl ^ci-ivip - g), P - q; 
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Notably, a strong momentum dependence of the integral 
kernel in Eq.(15) rules out the customary BGS-type so- 
lution A{p) = const. Goncomitantly, because of the van- 
ishing DOS of the Dirac liquid the Goulomb coupling 
would have to be strong enough to provide for a possible 
onset of pairing. 

At T = and with the velocity renormalization (13) 
taken into account, one can cast Eq.(15) in the following 
approximate form 



A(;,) = i> ^ 



q^dq 



|p-q|(A'7-F(7rgiV/8V2)|p-qh) 



A(9) 



(16) 

where 4> is the angle between the vectors p and q, and 
the additional factor of y/2 is characteristic of the typical 
energy- and momentum transfers [uj « vq) contributing 
to the integral. 

It is worth noting that, contrary to the naive expec- 
tation of its uniformly weakening effect on the effec- 
tive momentum-dependent coupling constant, the veloc- 
ity renormalization (13) affects the gap equation in a 
more subtle way. Namely, the above expectation is only 
valid for weak couplings, whereas in the strong coupling 
regime 17 3> 1 one finds a nearly complete cancellation 
between the velocity-dependent factors in the numerator 
and denominator of Eq.(16). 

In what follows, we analyze the gap equation along the 
lines of the earlier studies of chiral symmetry breaking—. 
In doing so, we primarily focus on its starkly novel fea- 
tures that are robust enough to sustain the approxima- 
tions involved. Namely, by following the same procedure 
as in Refs4i2iii^ we, first, linearize Eq.(16) with respect 
to the unknown mass function and then transform it into 
the differential form 



(fA{p) 
dp^ 
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where Srj = rj — rij\j and the effective coupling 
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interpolates between the weak {g w g) and strong 
{g sa 8^/2 /ttN) coupling limits, which are attainable by 
putting r/N — 0,g ^ and ijn = r/,g ^ 00, respectively. 
In order for Eq.(17) to be consistent with the origi- 
nal integral Eq.(15), it has be complemented with the 
boundary conditions 
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\p=A/v = 



[il + r,N)^{p)+P 



dp 



lp=A 



(19) 



The linearized Eq.(17) can also be interpreted as the 
radial Schroedinger equation for the zero-energy state 
in the two-dimensional spherically-symmetrical poten- 
tial U{r) (X l/r^~^^. This analogy suggests that for 
< 5r] < 1 the zero-energy solution might exist for all 
couplings g in excess of a finite threshold value Ijc- 

In order to illustrate the behavior of the physically rel- 
evant solution of Eq.(17) we employ the WKB method 
and consider a linear combination of two linearly inde- 
pendent functions 



A±(p) 



C± 



pl-5,,/2p(p)l/2 



exp(±i / P{p')dp') (20) 



where, for p > k = [(1 4- riN)/^^)]^^^''^^, 
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On the general grounds, one expects the physically ac- 
ceptable solution to monotonically increase with decreas- 
ing momentum, eventually leveling off at a finite value 
A = A(p)|p=o- Since the linearized equation (17) ceases 
to be applicable at momenta below p ^ A/w, a reliable 
use of Eqs. (20,21) would be limited to the range of pa- 
rameters p ^ A/v > K. 

Strictly speaking, the applicability of the WKB tech- 
nique requires the (quantized in multiples of tt) values 
of the momentum integral in the exponential of Eq.(20) 
to be large, in which case the mass function A{p) de- 
velops unphysical oscillations as a function of momen- 
tum. Therefore, in order to find a monotonic solution 
which satisfies the boundary conditions (19), we choose 
C+/C^ — i and equate the aforementioned integral to 
its lowest possible 'quantized' value 
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where 



^V5^b^''--^>) + ^A + <5A = .r (22) 
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1 + W -'5^/2(1- (.9c/.9)) 



8. = tan^'^^^^^^P^W^^. (23) 

At 7/ := Eq.(17) becomes scale-invariant and Eq.(20) 
amounts to the solution obtained in Refsi^i^ 
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A 



^^^^-^sin(-y2^^1n- + 2<5) (24) 



where 5 = (5a = <5a = tan~^ ^2g — 1. It is worth not- 
ing that, regardless of the formal restrictions on the ap- 
plicability of the WKB approximation, Eq.(24) remains 
accurate in the entire range A%p < A^. 

The maximum value of the mass function A = A(p = 
A/v) is given by Eq.(22) 



A = .Aexp(-^===) 



(25) 



Obviously, such a solution requires the coupling to be in 
excess of the critical value 



5c 



(26) 



For A^ = 4, Eqs. (18) and (26) yield the critical couphng 
gc = 1.13 which appears to be tantalizingly close to (al- 
beit still higher than) its actual value in graphene on 
the Si02 substrate, but below that in free-standing (sus- 
pended) graphene. Conversely, in the asymptotic g — > oo 
regime the critical number of fermion species below which 
the mass gets generated is Nc = 16^/2/ n sa 7.18. 

For ?7 <C 1 the analysis of Eq.(22) reveals an increase 
in the critical value 

ffc-^(l + W + (3^<577)'/3 + ...), (27) 

showing that a nonzero mass can emerge even in the case 
of a momentum-dependent coupling. 

The gap equation (15) can also be used to estimate the 
mean-field critical temperature. For this approximate 
calculation we neglect the weak (at large N) renormal- 
ization (13) by putting 77 = 0. The critical temperature 
can then be evaluated by solving the equation 



/ — ./2(gtanh— - gc) ~ 27r 
4 pV 2T, ^^^ 



(28) 



where pq is defined as tanh(upo/2rc) — l/g. As a re- 
sult, Eq.(28) yields the relation between the gap and the 
(mean-field) critical temperature 



Tr 



|ln(l-g,/g)|' 



(29) 



which is markedly different from the universal ratio 
Tc/A = e'-^/tt characteristic of the standard BCS case. 

Despite the highly unconventional relation (29), both 
A and Tc vanish rapidly upon approaching the threshold 
coupling. Appyling the above rough estimates to un- 
doped free-standing graphene one obtains the gap and 
the critical temperature of order A ~ Tc ~ lOe^e^^ ~ 
AmeV which is comparable to the gap ~ lOmeV reported 
in Ref.^ (an order of magnitude larger gap found in epi- 
taxial graphenei is likely to be induced by the SiC sub- 
strate, though). 

A finite chemical potential /z represents, depending on 
its sign, a finite density of either particles or holes. It 
gives rise to a mismatch between the occupation num- 
bers of the particles and holes which sets a lower limit of 
the momentum integration in Eq.(15). A straightforward 
analysis then shows that at p ^ fj,/v the mass function 
A{p) tends to level off at a value below that at /i = 0, 
whereas for p > /i/w the solution obtained for 11 — 
remains essentially intact. 

It also implies that the excitonic pairing should be sen- 
sitive to such particle-hole symmetry breaking factors as 
electron-hole separation ('puddles'). As follows from the 



above estimates, the presence of puddles with density 
n ~ 10^^cm~^ could already hinder the excitonic gap 
A ~ lOmeV. 

Electron-phonon interactions and Cooper pair- 
ing 

Next, we discuss putative Cooper pairing instabilites in 
graphene. For interactions dominated by small momen- 
tum transfers, the onset of the Cooper instability would 
be contingent with the existence of a range of parameters 
where the effective electron interaction becomes attrac- 
tive due to a more efficient screening of the Coulomb 
interaction. This situation can potentially occur, e.g., 
for a large-e substrate (e.g., if /O2) or in the presence of 
water molecules trapped between the graphene sheet and 
the substrate. 

The Cooper pairing scenario would provide for an al- 
ternative explanation of the data of Rcf.^^ indicating the 
presence of a ~ lOOmeV gap at the Fermi energy in epi- 
taxially grown graphene on a SiC substrate. The previ- 
ously proposed explanations of this finding include such 
mechanisms as electronic states in the underlying buffer 
layer, charging/band bending effects, coupling to the sur- 
face excitations, as well as a trivial finite size effect. As 
yet another idea, a phonon-mediated inelastic tunneling 
was argued to be behind the apparent TDOS suppression 
which was then misinterpreted as a gap openingi^. 

However, considering that the system appears to be 
naturally doped, /i ~ OAeV^, it is conceivable that at 
such a sizable doping the Coulomb repulsion would be 
rather strongly screened out and the electron-phonon 
coupling could then get a chance to become dominant, 
resulting in the Cooper instability. 

The phonon spectrum of graphene consists of a total 
of 6 modes, which list includes (in- and out- of phase) 
in-plane longitudinal/transverse and out-of-plane ones. 
Dispersion relations of the pertinent low-energy acoustic 
modes with momenta near the F point can be obtained 
from the (somewhat simplified for the sake of our quali- 
tative discussion) elastic energy of the graphene sheet 



F = ^muf + {d,hf 



-c^{diUj + djUi 



d.hdjhf] 
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which governs the in-plane phonons with a linear disper- 
sion {yii{q) = cq where c ^ lO^rn/s) and linear coupling 
to the Dirac fermions as well as the out-of-plane ('fluxor') 
modes with a quadratic dispersion (5^2(9) = t^q^) and 
quadratic couplingi^. The optical modes near the F- 
point, as well as all the modes at the i^i^-ones are 
gapped and their dispersions remain nearly constant 
(170(9) = f^o) at small deviations from the corresponding 
points in the momentum space. 

The corresponding matrix elements controling the 
electron-phonon coupling vary with momentum as 



M," ex ( 



7^0(9)' 



(31) 



where a = 0,1,2 and mc — 2 y. lQ~^^g is the mass of 
carbon atom. 

The overall effective electron interaction mediated by 
the phonon modes is given by the expression 

VpH{uj,ci)= J2 DaiLu,q)\M^\^ + 

a=QS 

Y.j%\K\"D2{^'+^c(^)D2{^'-^c{_) (32) 



where ijj'^ ~ uj' ± a;/2, q'_j_ = q' ± q/2, and 

^a{q) 



Daiuj,q) 



uj^^nl{q)+iO 



(33) 



is the propagator of the a-type mode. 

A relative importance of various phonon-induced elec- 
tron interactions can be readily assessed in the quasi- 
static (lo — 0) approximation. When treated this way, all 
the gapped and linear in-plane modes would have given 
risen to an effective interaction with the Fourier trans- 
form 



.(o.l),„^_ ^\M^\' 



C^'('Z) = -E 



a=0,l "^-^^ 



const 



(34) 
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which is short-ranged in the real space {V\' (r) ex 
— (5(r)). By contrast, the quadratic fluxor mode would 
produce a long-ranged potential with the Fourier trans- 
form 



K 



(2)^„^_ 



ph 



(^) = -E 



|M||2 



n2ik + q) + n2{k) 



r{2). 
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(35) 



which decays with the distance as V\ (r) ex — l/r . 

Thus, depending on the relative strength of this cou- 
pling, the fluxor mode could provide for the strongest 
tendency towards the Cooper pairing, as compared to 
the gapped and linear acoustic modes. 

Focusing solely on the electron-phonon coupling, one 
can derive the Cooper analog of the gap equation (15) for 
the 'on-sheir self-energy Asc'ip) = ^a{vp ~ /i, p). Upon 
integrating over the frequencies one obtains 



^sc 



b)-E 



|Mp-qP 



Qaip - cO + Eiq) 2E{q) 



^-('^)tanh« 

2r 



(36) 
Contrary to the case of the excitonic instability at /x = 0, 
the solution of Eq.(36) exists for an arbitrarily weak 
electron-phonon interaction (which, however, is assumed 
to be stronger than the screened Coulomb one). This so- 
lution appears to be essentially independent of the mo- 
mentum for p ^ ^/v, and the gap conforms, by and large, 
to the standard BCS formula 



^sc« E 

a=0,l,2 



Eae 



-1/A„ 



(37) 



where the characteristic energy scales and dimensionless 
couphng constants corresponding to the different modes 
are given by the expressions 



Eq ^ minlflo, fi]. 
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and 
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respectively. Here D « leV/A is the deformation poten- 
tial, a — \A2A is the lattice spacing, and JIq ~ 0.2eV is 
the energy of the Raman-active gapped phonon modes. 

Contrary to excitonic pairing, the Cooper one is facili- 
tated by an increasing chemical potential (of either sign) , 
while it gets suppressed by a Zeeman-coupled magnetic 
field. Also, the logarithmic factor in the coupling to the 
fiuxor mode A2 is another manusfestation of the long- 
ranged nature of the effective electron interaction (35) 
mediated by this mode. 

As far as the actual values of the couplings (39) are 
concerned, those predicted by the density functional 
theoryi^ seem to be much too weak to explain the ex- 
perimentally observed downward renormalization of the 
Fermi velocityis^ii^. Specifically, apart from being of 
the sign opposite to that of the predicted Coulomb-only 
renormalization (13), the experimentally measured ve- 
locity change appears to be quite large in magnitude 
(« 20%) and suggestive of the effective elcctron-phonon 
couphng of order Aq ~ 0.2 — 0.4. 

The analysis of Ref,^^ limited to the gapped in-plane 
modes indicates that, depending on the mode's symme- 
try, the corresponding coupling constant might become 
strongly enhanced at low momenta by the Coulomb inter- 
actions. Namely, in Ref.^^ the coupling constant for the 
optical 'Al'-mode at the iiT-point was found to increase 
by over an order of magnitude from its microscopic value 
as a result of the scaling Ao(p) — \o{v{p)/v)'^. 

For the largest estimated values of the above parame- 
ters the gap can be as high as A w OAeVe'^'^ ^ 30meV 
which compares favorably with the experimentally ob- 
served ~ lOOmeV gapi^. 

The lack of a significant momentum dependence in 
Eq.(36) would give rise to the BCS relation between Asc 
and the mean- field critical temperature T^'~^ . Although 
the above estimates suggest that for the experimentally 
relevant electron densities n ^ 10^^cm~^ the latter could 
be as high as T^'^ ~ 150K, the prospects of observing the 
Cooper instability in graphene are likely to be hampered 
by a number of factors. 

Firstly, the above Peierls- and Cooper-pairings would 
primarily manifest themselves through a "pseudogap" 
in the electron spectrum which emerges below the cor- 
responding mean-field critical temperatures. However, 
a true transition in such a two-dimensional system as 



graphene would only take place at a temperature deter- 
mined by the (un)binding of the vortex-antivortex exci- 
tations of the order parameter in question. 

The critical temperature of the underlying Kosterlitz- 
Thouless (KT) transition is related to the rigidity with 
respect to the fluctuations of that order parameter. The 
latter can be computed by expanding the pertinent sus- 
ceptibility 



n.j(o,g ), 



^'E 



duj 



1 9^0 
1 



^^TrT,Giu;,p + q)TjG{u;,p) « ^A (40) 



to second order in the center-of-mass momentum q, which 
procedure yields the KT temperature 



T, 



KT 



^Ps 



I- 



(41) 



The observation that the KT transition temperature is 
smaller than the mean-field one by an order of magnitude 
is in agreement with the previous work of Refsi^. 

Albeit being robust against potential disorder, the on- 
set of the Cooper pairing can be hindered by a rippling 
of the graphene sheet. The latter was argued to pro- 
duce effective magnetic fields as high as B ^ 1 — ST^i, 
whose pair-breaking effect could potentially eliminate 
the propensity towards the Cooper instability altogether, 
since for Asc ^ O.leV^ the upper critical field Bc2 — 
^oAgfj/2TTv'^ would be of order -- 4T. 

Conclusions 

In summary, we carried out the analysis of excitonic 
(Peierls) and Cooper pairing instabilities in graphene 
which presents a physically relevant example of the nodal 
fermion systems with a power-law DOS and long-range 
pairwise interactions. 

We find that undoped graphene can be prone to the on- 
set of a variety of excitonic instabilities, provided that the 
Coulomb coupling is sufficiently (albeit not unphysically) 
strong. Unlike in the low-energy (hence, weak-coupling) 
BCS-like pairing scenaria of Refs.^'"^, the non-BCS solu- 
tion of the gap equation Eq.(15) is strongly momentum- 
dependent and requires the coupling to be in excess of 
a finite threshold value gc- Moreover, the perturbative 
downward renormalization of the running coupling con- 
stant g{p) does not appear to have a dramatic impact on 
the onset of this excitonic instability (or a lack thereof). 

With increasing doping, the excitonic instability be- 
comes suppressed, while a propensity towards the zero- 
threshold Cooper pairing driven by the phonon-mediated 
attraction gets stronger (in practice, however, it remains 
countered by the residual screened Coulomb repulsion). 
Should either of these instabilities fully develop, it would 
give rise to a momentum-dependent gap function and/or 
other experimentally testable non-BCS features. 

A better understanding of the phenomena resulting in 
a finite gap in the nodal fermion spectrum is going to be 
instrumental for the proposed applications of graphene 



in post-silicon electronics. Therefore, one should expect 
that future experiments will help to ascertain the actual 
status of the above scenarios before long. 
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This research was supported by NSF under Grant DMR- 
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Note added: 
After this work was first made public (|arXiv:0807.0676p . 
the results of two new Monte Carlo simulations were 
released, according to which the aforementioned exci- 



tonic transition in free-standing graphene might occur 
at the critical coupling as low as ^c = 1-1 for iV = 43^ 
and the number of fermions as high as Nc — 9.4 for 
g — > ooS2, (note that in both preprints the quoted num- 
ber of fermions N — 2 pertains to the /owr-coniponent 
bi-spinors, which corresponds to iV = 4 in this paper). 
Also, in the recent work^** the integral gap equation (10) 
was solved numerically without any further approxima- 
tions and the critical number of fermions was found to 
be Nf, « 7. All the above values appear to be remarkably 
close to those obtained in this paper. 
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